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ABSTRACT 

We report the discovery and monitoring of radio emission from the Type Ic SN 2002ap ranging in frequency 
from 1.43 to 22.5 GHz, and in time from 4 to 50 days after the SN explosion. As in most other radio SNe, the 
radio spectrum of SN 2002ap shows evidence for absorption at low frequencies, usually attributed to synchrotron 
self-absorption or free-free absorption. While it is difficult to discriminate between these two processes based on 
a goodness-of-fit, the unabsorbed emission in the free-free model requires an unreasonably large ejecta energy. 
Therefore, on physical grounds we favor the synchrotron self-absorption (SSA) model. In the SSA framework, at 
about day 2, the shock speed is « 0.3c, the energy in relativistic electrons and magnetic fields is ~ 1.5 X 10 45 erg 
and the inferred progenitor mass loss rate is « 5 x 10~ 7 Mq/vt (assuming a 10 3 km sec" 1 wind). These properties 
are consistent with a model in which the outer, high velocity supernova ejecta interact with the progenitor wind. 
The amount of relativistic ejecta in this model is small, so that the presence of broad lines in the spectrum of a 
Type Ib/c supernova, as observed in SN 2002ap, is not a reliable indicator of a connection to relativistic ejecta and 
hence 7-ray emission. 

Subject headings: radio continuum:supernovae — supernovae:individual(SN2002ap) — starsxircumstellar 
matter — gamma rays:bursts 



1. INTRODUCTION 

Type Ib/c supernovae (SNe) enjoyed a broadening in interest 
over the last few years since their compact progenitors (Helium 
or Carbon stars) are ideal for detecting the signatures of a cen- 
tral engine. Such an engine is expected in the collapsar model 
(Woosley 1993; MacFadyen, Woosley & Heger 2001), the cur- 
rently popular model for long-duration 7-ray bursts (GRBs). In 
this model, the engine (a rotating and accreting black hole) pro- 
vides the dominant source of explosive power. The absence of 
an extensive Hydrogen envelope in the progenitor star may al- 
low the jets from the central engine to propagate to the surface 
and subsequently power bursts of 7-rays. 

Separately, the Type Ic SN 1998bw (Galama et al. 1998) 
found in the localization region of GRB 980425 (Pian et 
al. 2000) ignited interest in "hypernovae" 1 . Regardless 
of the controversy over the association of SN 1998bw with 
GRB 980425, or equivalently, the controversy over the relation 
of the extremely underluminous GRB 980425 to the cosmolog- 
ical GRBs, one fact is not in dispute: SN 1998bw is a most in- 
teresting SN. First, the SN exhibited the broadest absorption 
lines to date, about 60,000 km sec" 1 (Iwamoto et al. 1998; 
Woosley, Eastman & Schmidt 1999). Second, modeling of the 
optical spectra and lightcurves suggested a large energy release, 
£51 r-j 30, where E$\ is the SN energy release in units of 10 51 
erg. Finally, and most relevant to the issue of GRB connec- 
tion, the SN was the brightest radio SN at early times; robust 
equipartition arguments led to an inferred energy of £p J> 10 49 
erg in ejecta with relativistic velocities, T ~ few (Kulkarni et al. 



1998, hereafter K98). Until SN 1998bw, no other SN showed 
hints of such an abundance of relativistic ejecta. Tan, Matzner 
& McKee (2001) explain the relativistic ejecta as resulting from 
an energetic shock as it speeds up the steep density gradient of 
the progenitor. The 7-ray and radio emission would then arise 
in the forward shock. 

From the perspective of a GRB-SN connection, what matters 
most is the presence of relativistic ejecta. 7-ray emission traces 
ultra-relativistic ejecta, but as was dramatically demonstrated 
by SN 1998bw, the radio serves as an equally good proxy for 
relativistic ejecta with the added advantage that the emission 
is not beamed. Given this, we began a systematic program of 
investigating at radio wavelengths all Ib/c SNe with features 
similar to SN 1998bw: a hypernova or broad optical lines. 

Y. Hirose discovered SN2002ap in M74 (distance, d ~ 
7.3Mpc; Smartt et al. 2002) on 2002, Jan. 29.40 UT (see 
Nakano 2002). Mazzali et al. (2002) inferred an explosion date 
of 2002, Jan. 28 ± 0.5 UT. Attracted by the broad spectral fea- 
tures (e.g. Kinugasa et al. 2002; Meikle et al. 2002) we began 
observing the SN at the Very Large Array (VLA 2 ). 

2. OBSERVATIONS 

We observed SN 2002ap starting on 2002, February 1 .03 UT. 
and detected a radio source coincident with the optical position 
at a(J2000)= 01 h 36 m 23.92 s , 5(J2000)=+15°45'12.87", with a 
l-cr uncertainty of 0.1 arcsec in each coordinate (Berger et al. 
2002a). A log of the monitoring observations and the resulting 
lightcurves can be found in Tab. |l]and Fig. [I], respectively. 



'There is no accepted definition for a hypernova. Here we use the term to mean a supernova with an explosion energy significantly larger than 10 5 ' erg. 
2 The VLA is operated by the National Radio Astronomy Observatory, a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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2. 1 . The Radio Spectrum of SN 2002ap 



The peak radio luminosity of SN2002ap, L p (5GHz) - 3 x 
10 25 erg sec -1 Hz" 1 , is a factor of 20 lower than the typical 
Ib/c SNe (Weiler et al. 1998), and - 3 x 10 3 times lower than 
SN 1998bw (K98). The time at which the radio emission peaks 
at 5 GHz is t p ~ 3 day, which may be compared with 10 days 
for SN 1998bw (K98), and 10-30 days for the typical Ib/c SNe 
(Weiler et al. 1998; Chevalier 1998, hereafter C98). 

The spectral index between 1.43 and 4.86 GHz, /3?jf, 
changes from ~ 0.5 before day 6, to ~ -0.3 at t sa 15 days 
(F„ oc v@), while holds steady at a value of w -0.9. This 
indicates that the spectral peak, v p , is initially located between 
1.43 and 4.86 GHz, and decreases with time. This peak could 
be due to synchrotron self-absorption (SSA) or (predominantly) 
free-free absorption (FFA) arising in the circumstellar medium 
(CSM). Regardless of the dominant source of opacity, the emis- 
sion for frequencies v > v p is from optically-thin synchrotron 
emission. 

Massive stars lose matter via strong stellar winds through- 
out their life and as a result their CSM is inhomogeneous with 
density, p(r) oc M w v~ l r~ 2 . Here, r is the distance from the star, 
M is the rate of mass loss, and v„ is the wind speed, which is 
comparable to the escape velocity from the star. The progeni- 
tors of Type II SNe are giant stars which have low v„ ~ 10 km 
s . Consequently the CSM is dense and this explains why the 
FFA model has provided good fits to Type II SNe (e.g. Weiler, 
Panagia & Montes 2001). 

On the other hand, the progenitors of Type lb and Ic SNe 
are compact Helium and Carbon stars which have high escape 
velocities and therefore fast winds ~ 10 3 km sec -1 . Thus, a pri- 
ori, the CSM density is not expected to be high. C98 reviews 
the modeling of radio lightcurves of Ib/c SNe and concludes 
that there is little need to invoke free-free absorption. However, 
synchrotron self absorption is an inescapable source of opacity 
and must be included in the modeling of Type Ib/c SNe (C98; 
K98). 

Low frequency observations provide the simplest way to dis- 
criminate between the two models. In the SSA model, the peak 
frequency is identified with the synchrotron-self absorption fre- 
quency, v a , and F v (v <; u a ) oc v 5 l 2 . In the FFA model, the 
free-free optical depth is unity at the peak frequency, vjf and 
F v (v < i/ff) decreases exponentially. Lacking the requisite dis- 
criminatory low frequency data we consider both models. 

2.2. Robust Constraints 

Before performing a detailed analysis, we derive some gen- 
eral constraints using the well-established equipartition argu- 
ments (Readhead 1994; K98). The energy of a synchrotron 
source with flux density, F p {v p ), can be expressed in terms 
of the equipartition energy density, U /U eq = +rj~ ll )/2, 
where rj = 6> s /6> eq , 6 eq « UOd^ 1 F^ y z£g^ 35)/34 Mas, and 

U eq = LI X 10 ^Mpc^p.mJy^p.GHz^eq./jas er S- 

At about day 7, F p (v p = 1.4GHz) w 0.3 mJy (see Fig. |l]). 
Thus, 9 eq (t = 7d) w 40 /ias, or r w 4.5 x 10 15 cm. The result- 
ing equipartition energy is E eq i=s 10 45 erg, the magnetic field is 
B eq ~ 0.2 G, and the average velocity of the ejecta is v e q ~ 0.3c 
as inferred above. We note that any other source of opacity (e.g. 
free-free absorption) would serve to increase 8 eq , E eq , and v eq . 



3. A SYNCHROTRON SELF-ABSORPTION MODEL 

The synchrotron spectrum from a source with power-law 
electron distribution, N(j) oc r y~ p for 7 > j m \ n (§2. 1 ) is 



F v =F vfi (vv Q ) 1 (l-e »)— -— -, (1) 

Fi{vo,v m ,p) F 2 (v,v mi p) 

where the optical depth at frequency v is given by 



T „ = T -o(^o)- (2+ " /2) 



F^jv.VnuP) 
F2(fO,^m,P)' 



and 



Ft{v,v m ,p)= / F{x)x {p - l) l 2 dx; 



(2) 



(3) 



see Li & Chevalier (1999; hereafter LC99). Here x„, = vjv m 
(see Rybicki & Lightman 1979), and v m is the characteristic 
synchrotron frequency of electrons with 7 = 7 m ; n . The subscript 
zero indicates quantities at a reference frequency which we set 
to 1 GHz. Finally, v a is defined by the equation r Ua = \. 

The evolution of the synchrotron emission depends on a 
number of parameters. Following C98, we assume that p, and 
the fraction of energy in electrons (e e ) and magnetic fields (e#) 
in the post-shock region remain constant with time; unless oth- 
erwise stated, e e = eg = 0.1. The evolution of the synchrotron 
spectrum is sensitive to the expansion radius of the forward 
shock front, r s oc t m , which is related to the density structure 
of the shocked ejecta and that of the CSM. We allow for these 
hydrodynamic uncertainties by letting F Vj q oc t^ r and tq oc t? T , 
where is the time days since the SN explosion. In the model 
adopted here, both these indices depend on m and p. It can 
be shown that the temporal index of the optically thin flux, 
a = ctF + ol t . The synchrotron characteristic frequency, v m , 
is particularly useful for inferring the CSM density, and we 
parametrize it as v m = v m ,ot^ v '" GHz where v m $ = v m GHz. For 



typical values of m, and p(r) otr ,a v 



-0.9. 



With these scalings and Eqs. [l]-|3| we carry a least-squares fit 
to the data. Given the lack of early optically-thick data (i.e. 1 .43 
GHz) it is not surprising that our least-squares analysis allows 
a broad range of values for a T . In Fig. [j] we plot fits span- 
ning the minimum x 2 '- a T = —1.3,— 2.1,— 3 (corresponding to 
X 2 = 40,43,46, respectively and 21 degrees of freedom). We 
note that for other Ib/c SNe a T range from -2 to -3 (C98; 
LC99). 

The fits, in conjunction with Eqs. 13-15 of LC99 allow us 
to trace the evolution of the size of the expanding blastwave 
(r s ), the total (magnetic+electrons) energy (E), and the electron 
density (n e ) in the shock (Fig. ||). We find that for a T = -1.3, 
r s oc f 0,25 i.e. the blastwave decelerates. However, a T = -3 
provides the expected r s oc t . Adopting this physically rea- 
sonable model, we obtain: Tnit) = 1.2 x 10 3 f/, F ufi (t) = 2.9fp 
/zJy, and p = 2. From Fig. || we note that the early shock ve- 
locity is high, 0.3c, regardless of the choice of a T , and clos e to 
that derived from the simple equipartition arguments (§2.2). 

A simple extrapolation of our model, with a cooling break 
between the radio and X-ray bands (as expected from the 
value of B eq inferred in §2.2) yields an X-ray flux compara- 
ble to the measurement by XMM-Newton on 2002, Feb 3, 
Fx ~ (1.2 ±0.3) x 10~ 14 erg cm" 2 sec" 1 (Rodriguez Pascual 
et al. 2002; Soria & Kong 2002). In addition, the hardness 
ratio of the X-ray emission is consistent at the 2a level with 
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p = 2 inferred from the radio data. Thus X-ray emission due 
to other mechanisms (e.g. thermal brehmsstrahlung) must be 
quite small. 

The mass loss rate of the progenitor star is estimated from 
r s and n e , M w = 8Tr£n e m p r 2 v w « 9 x 10~ 9 ^~°q 8 M q yr" 1 , where 
the compression factor is £ = 1 /4, the nucleon-to-electron ratio 
is taken to be 2 and v w = 10 3 km sec -1 . Knowing B eq and our 
assumed e e we find v m ~ 10 7 Hz and thus M w w5 x 10~ 7 Mq 
yr" 1 - similar to that inferred for SN 1998bw (LC99). 

There are two consistency checks. First, with this M w , free- 
free absorption is negligible. Second, the kinetic energy of the 
swept-up material is 2 x 10 46 erg - consistent with our estimate 
of the equipartition energy and e e . 

3.1. The SSA Model in the Context of a Hydrodynamic Model 

The results of §|| can be tied in to a fairly simple hydro- 
dynamic model. Matzner & McKee (1999) show that for 
the progenitors of Vole SNe (compact stars with radiative en- 
velopes) the ejecta post-explosion density profile can be de- 
scribed by power laws at low and high velocities, separated 
by a break velocity, v ej , b = 5150(£ 5 i/Mi) 1/2 « 2 x 10 4 km 
sec -1 ; where the mass of the ejecta is M e j = \QM\Mq. We 
use E 51 « 4-10 and M\ « 0.25-0.5 for SN2002ap (Maz- 
zali et al. 2002). At v s w 0.3c, the density profile is given 
by p s « 3 x 10 % £ 3 i 59 M7 259 r 3 v- 1018 g cm" 3 . This profile ex- 
tends until radiative losses become important when the shock 
front breaks out of the star. Using Eq. 32 of Matzner & Mc- 
Kee (1999) this happens for v s « 1 .5c (assuming a typical 1 R 
radius for the progenitor star). Thus, the outflow can become 
relativistic. 

Using the self-similar solution of Chevalier (1982) the ve- 
locity of the outer shock radius, R, (assuming a p = Ar~ 2 CSM) 
is 

- = 52, 300£ 5 Y 4 M7°- 32 A; 01 V 12 km sec" 1 , (4) 

where A, = (M.,/10 -5 M yr 1 )(v vv /10 3 km sec" 1 )" 1 . The 
shock velocity, R is insensitive to the circumstellar wind den- 
sity. Thus, we find that the velocities inferred from the radio 
observations of SN 2002ap can be naturally accounted for by 
the outer supernova ejecta. 

The energy above some velocity V is 

r°° l 

E(v > V) w / -p f v 2 4irv 2 t 3 dv = 7.2 x lO 44 ^ 59 ^ 2 ' 59 ^ 5 ' 18 ergs, 
J v 2 

(5) 

where V5 is the velocity in units of 10 5 km s 1 . For the pre- 
ferred SN 2002ap parameters, E(v > V) « 4.2 x lO 48 ^ 518 erg. 
There is therefore plenty of energy in the high velocity ejecta to 
account for the observed radio emission. 

Indeed, given the over-abundance of E(v > V5) relative to the 
energy inferred from the radio emission, we wonder how secure 
the are the estimates of £51 and Mi of Mazzali et al. (2002). 
In particular, £51 and Mi are derived from early time optical 
observations and are certainly subject to asymmetries in the ex- 
plosion. For SN 1998bw, the asymmetric model of Hoflich, 
Wheeler & Wang (1999) yielded £51 ~ 2, an order of magni- 
tude smaller than that obtained from symmetrical models (e.g. 
Iwamotoetal. 1998). 

3.2. Interstellar Scattering & Scintillation 



Interstellar scattering and scintillation (ISS) is expected for 
radio SNe (see K98). Indeed, the perceptible random devia- 
tions from the model curves (see Fig. [j]), which account for the 
high Xmin cou ld arise from ISS. 

We compute the ISS variability, specifically the modulation 
index (the ratio of the rms to the mean) of the observed flux, 
using the ISS model of Goodman (1997). Along the direction 
to SN 2002ap using the Galactic free electron model of Taylor 
& Cordes (1993) we find a scattering measure, SM_3 5 sa 0.65 
m -20 ' 3 kpc and the distance to the scattering screen of d sa - w 0.5 
kpc. 

For radio SNe there are two basic regimes of interest sep- 
arated by the so-called transition frequency, v a . For v > v tI , 
the weak scattering regime, the modulation index is „ = 
(t'tr/V) 17 / 12 [for 9 S = r s /d smaller than the Fresnel size, 6-p = 
8.1(c/ scr i/GHz) -1 ^ 2 /^as] and m UA = m u ^{9'p/0 f ) 1 l ( ' otherwise. 
For v < v tI , the strong (refractive) scattering regime, m„R = 

(!//^tr) 17/30 [for 6 S < 0r = 6> F , tl -(t'tr/^) 11/5 , where <V is the Fres- 
nel size at v tI ], and m„, q = w^r^r/^) 7 ' 6 otherwise. 

For SN 2002ap we find v a « 7.3 GHz. From Fig. || the source 
angular diameter, over the period 5-20 days, is 9 S « 60- 160 
/ias, from which we estimate the following values: m%M ~ 5%, 
7«4.86 ~ 10%, and W1.43 w 40%. 

We estimate the actual modulation index empirically by 
adding m v F v in quadrature to each measurement error so that 
the reduced x^m * s ur, ity- Here F u is the model flux described 
in §|[ We find mg.46 ks 0.1, m^gg ~ 0.2, and m\ .43 w 0.3, in 
good agreement with the theoretical estimates. This provides 
an independent confirmation of the size, and hence expansion 
velocity of the ejecta. 

4. A FREE-FREE ABSORPTION MODEL 

In this model, the spectrum is parametrized as (Chevalier 
1984; Weiler et al. 1986): 

F v =K x vft P d e-^ 
T v =K 2 v? l t 5 d , (6) 

where 1/5 = 5v GHz. We find an acceptable fit (\ 2 = 40 for 21 
degrees of freedom) yielding: K\ m 2 mJy, K2 ~ 0.4, a w -0.9, 
j3 w -0.9, and 5 w -0.8. With these parameters and Eq. 16 of 
Weiler et al. 1986) we find M w w 5 x 10" 5 M Q yr" 1 for v w = 10 3 
km sec" 1 . 

In this model, one day after the explosion ss 3.2 GHz, 
and /v(t'ff) « 1.1 mJy (Fig. [j]). The unabsorbed flux at the 
peak of the synchrotron spectrum is F v (v a ) w 3(z^/3.2GHz)" ' 9 
mJy (note v a < u/f in the FFA model) for which r eq « 7.5 x 
10 15 (j/ fl /3.2GHz)- 3/2 cm. Thus v eq w 3c(iy fl /3.2GHz)" 3 / 2 , 
which corresponds to T = 2(z/ fl /3.2GHz) _I if relativistic effects 
are taken into account (R. Sari priv. comm.). Alternatively, if 
we fix the expansion velocity to the optical value, v, w3x 10 4 
km sec" 1 (Mazzali et al. 2002), we find a brightness temper- 
ature, Tb « 4 X 10 13 K — clearly in excess of the equiparti- 
tion temperature, again necessitating a high bulk Lorentz factor, 

r ~ 10 2 . 

Thus, even if v a = v<s (in which case free-free opacity would 
not be necessary in the first place), the FFA model requires 
truly relativistic ejecta, or alternatively a large departure from 
equipartition, resulting in E w 7 x 10 50 (^ fl /3.2GHz)~ 9 erg (for 
v s w 0.5c instead of 3c). Clearly, the energy requirement would 
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increase by many orders of magnitude if v a <C vr. 

5. DISCUSSION AND CONCLUSIONS 

The type Ic SN 1998bw, likely associated with GRB 980425, 
was peculiar in two ways. It exhibited broad photospheric ab- 
sorption lines, and exceedingly strong radio emission at early 
times. The latter was interpreted to arise from relativistic (V ~ 
few) ejecta of total energy Er > 10 49 erg. These two pecu- 
liarities made sense in that the simple theory suggested that 
broad photospheric features are a reliable indicator of relativis- 
tic ejecta, a necessary condition for 7-ray emission. 

The type Ic SN 2002ap elicited much interest because it too 
displayed similar broad lines. However, from our radio obser- 
vations we estimate the energy in relativistic electrons and mag- 
netic fields to be quite modest: E w 1.5 x 10 45 ergs in ejecta 
with a velocity w 0.3c. In view of this, the absence of 7-rays 
from SN2002ap is not surprising (Hurley et al. 2002). Both 
the energy and speed of the ejecta can be accounted for in the 
standard hydrodynamical model. We thus conclude that broad 
photospheric lines are not good predictors of relativistic ejecta. 

Separately, the broad photospheric features led modelers to 
conclude that SN 2002ap was a hypernova with an explosion 
energy of £51 ^4-10 erg (Mazzali et al. 2002). However, 
the radio observations paint a different picture. The low Er 



make SN 2002ap an ordinary Type Ib/c SN (or perhaps even a 
low energy event; §2.1). We also draw attention to the signif- 
icant role played by asymmetries in forming the optical lines 
which, when not properly accounted for can lead to discrepant 
estimates (cf. the large variations in the estimated SN energy 
release for SN 1998bw; see Hoflich, Wheeler & Wang 1999; 
Iwamoto et al. 1998). Thus, we wonder how reliable are the 
inferred £51 and M\ values. Along the same vein, we note that 
Kawabata et al. (2002) suggest, based on spectro-polarimetric 
observations, a jet with a speed of 0.23c and carrying 2 x 10 51 
erg. Such a jet, regardless of geometry, would have produced 
copious radio emission. 

We end with three conclusions. First, at least from the per- 
spective of relativisic ejecta, SN2002ap was an ordinary Ib/c 
SN. Second, broad photospheric lines appear not to be a good 
proxy for either an hypernova origin or 7-ray emission. Third, 
radio observations offer a practical and accurate proxy for rela- 
tivistic ejecta. 
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Table 1 

Radio Observations of SN 2002ap 



Epoch 


fl.43±(T 


F4.86 ± °" 


^8.46 ± CT 


F22.5 ± 


(UT) 


OJy) 




(MJy) 


(A*Jy) 


2002 Feb 1.03 






374 ±29 




2002 Feb 1.93 


211 ±44 


384 ±50 


255 ±44 


348 ±165 


2002 Feb 2.79 


250 ±72 


453 ±50 


201 ±47 




2002 Feb 3.93 


410±41 


365 ±38 


282 ±34 




2002 Feb 5.96 


243 ±43 


262 ±48 


186 ±42 


170±91 


2002 Feb 8.00 


235 ±31 


282 ±32 


140 ±27 




2002 Feb 11.76 


337 ±68 


217±45 


1 1 1 ± 27 




2002 Feb 13.94 


292 ±38 








2002 Feb 18.95 


266 ±42 








2002 Mar 4.85+1 1.83 


157 ±34 








2002 Mar 18.77+19.97 


57 ±33 




25 ±25 





NOTE.— The columns are (left to right), (1) UTdate of each observation, and flux density andrms noise at (2) 1.43 GHz, (3)4.86 GHz, (4) 8.46 GHz, and (5) 22.5 
GHz. Observations with more than one date have been co-added to increase the signal-to-noise of the detection. 
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Time Since SN (days) 



Fig. 1. — Radio lightcurves of SN2002ap. The thick solid lines are our three synchrotron self-absorption models described in §[3[ 
with t v oc t~ l 3 , t v oc ?~ 2 1 , and t„ oc t " 3 in order of decreasing thickness. The dashed line is the model-fit based on free-free absorption 
(§Q). At 4.86, 8.46, and 22.5 GHz, the SSA and FFA models provide the same fit, since the opacity processes do not influence the 
optically-thin flux. The models diverge in the optically-thick regime, which underlines the importance of rapid, multi-frequency 
observations. 
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Fig. 2. — Inferred physical parameters based on the synchrotron self-absorption models described in §[| The panels are (a) time 
evolution of the total energy, (b) radius of the radio photosphere, (c) electron density in the shock as a function of radius, and (d) 



velocity of the shock front as a function of time. Models with t v oct , t v oc t , and r„ oc t are shown in order of decreasing 
thickness. The most likely fit is the one following r oc f° 9 (i.e. the model with t v oc f" 3 ). 



